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ABSTRACT 


The  following  Is  a  summary  of  research  activities  for  the 
period  October  1,  1975  to  September  30,  1980,  related  to  grant 
AF-AFOSR-76-2943.  The  primary  purpose  of  the  research  was  the  study 
of  forced,  thermoelastic  motion  of  plates  due  to  Laser  irradiation. 
This  summary  describes  the  theoretical  and  experimental  phases  of  the 


research  project. 
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1.0  Introduction 

One  of  the  primary  missions  of  the  U.S.  Air  Force  is  to  provide 
leadership  in  weapon  system  development,  and  to  be  cognizant  of  any  and 
all  advances  in  pure  and  applied  scientific  knowledge  which  has  a  bearing 
on  changes  in  the  state  of  the  art.  Current  developments  centering  on 
Laser  technology  have  advanced  to  a  stage  where  Lasers  must  be  considered 
as  potential  weapons.  A  high  power  Laser  has  the  ability  to  concentrate 
a  short  duration,  high  energy  flux  in  a  very  narrow  beam.  Thus  it  has  the 
ability  to  deposit  focused  radiant  energy  upon  opaque  solids  or  struc¬ 
tural  components.  Unlike  other,  conventional  weapons,  it  is  not  dependent 
upon  a  ballistic  path,  but  it  is  inherently  a  line  of  sight  weapon.  Because 
radiation  travels  with  the  speed  of  light  (approximately  3x10^®  cm/sec)  such 
a  weapon  can  inflict  almost  Instantaneous  damage  to  any  target  upon  which 
it  is  trained.  These  characteristics  obviate  the  need  for  a  ballistic 
computer  as  well  as  complex  lead  aiming  devices,  and  the  time  between 
target  acquisition  and  radiation  impact  is  negligible. 

The  present  study  was  concerned  with  the  interaction  of  a  laser  beam 
with  the  skin  of  an  aircraft,  a  re-entry  vehicle,  or  a  satellite,  i.e., 
a  thin  plate  or  shell.  Laser  interaction  with  plates  or  shells,  in  general 
is  a  complex  phenomenon.  For  purposes  of  experimentation  and  associated 
analysis,  it  is  possible  to  define  three  types  of  interaction  resulting 
in  (more  or  less)  separated  effects: 

(a)  Sudden  deposition  of  thermal  energy,  without  a  change  in  phase. 
This  causes  sudden  thermal  stresses  in  the  irradiated  plate.  Because  of  the 
rapidity  of  the  energy  deposition  process,  there  will  be  thermally  generated 


stress  waves 


(b)  Surface  vaporization  of  a  very  thin  layer  of  material,  plasma 
production,  plasma  heating,  shockwave  formation  in  plasmas,  etc.  This 
mode  of  interaction  results  in  suddenly  applied  surface  pressures  of  con¬ 
siderable  magnitude  to  the  solid  plate,  inducing  time-dependent  stress  waves 
and  deformations. 

(c)  Complete  local  vaporization  of  the  material,  and  the  resulting 
creation  of  openings  (holes  or  large  surface  cavities). 

It  is  to  be  noted  that  conditions  (a),  (b) ,  and  (c)  usually  coexist, 
but  they  can  be  (nearly)  separated  by  proper  choice  of  laser  and  target 
parameters. 

Laser  technology  and  laser  weapons  development  is  an  on-going  activity 
in  the  U.S.A.  and  in  other  countries.  By  virtue  of  its  mission,  the  U.S. 

Air  Force  must  keep  abreast  of  such  developments,  and  must  be  prepared  to 
evaluate  accurately  the  potential  and  destructive  capabilities  of  such 
weapons  when  and  if  they  become  operational.  In  addition,  detailed  know¬ 
ledge  of  the  capabilities  of  such  a  weapon  system  will  undoubtedly  suggest 
ways  and  means  to  either  avoid,  circumvent,  or  reduce  its  destructive 
effects  upon  potential  targets. 
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2.0  Sununary  of  Research  Accomplishments 
2.1  Theoretical 

Detailed  and  comprehensive  theoretical  studies  were  performed  to 
model  the  motion  of  a  rectangular  plate  when  subjected  to  laser  irradiation 
(see  references  2  and  4).  In  these  studies,  we  utilized  the  Bechtel  model 
for  Laser  heating. 

The  first  part  of  our  study  deals  with  the  thickness-stretch  motion 
of  a  transversely  constrained,  irradiated  slab.  This  part  of  the  study 
models  the  initial  response  of  the  plate.  The  initial  motion  of  the  plate 
is  predominantly  thickness-stretch  in  the  vicinity  of  the  irradiated  area. 
In  this  area,  stresses  and  displacements  are  primarily  in  the  direction  of 
plate  thickness.  However,  these  stresses  and  displacements  are  relatively 
small  compared  to  the  stresses  and  displacements  resulting  from  the 
gross  (predominantly  flexural)  motion  of  the  plate.  It  is  shox^n  that  in 
the  vicinity  of  the  Laser  beam,  a  dilatational  wave  is  set  up  which  moves 
in  the  direction  of  the  plate  thickness.  This  wave  causes  periodic  tension 
and  compression  stresses  on  the  median  surface  of  the  plate  (see  Fig.  1). 

The  second  part  of  this  study  is  concerned  with  the  gross-motion  of 
the  rectangular  plate  which  is  assumed  to  be  simply  supported  along  its 
boundaries.  The  plate  surface  is  irradiated  by  a  Laser  beam  at  an  arbi¬ 
trary  point.  Three  different  theories  are  used  to  model  the  time- 
dependent,  thermoelastic  motion  of  the  plate: 

(a)  Three-dimensional  elasticity  theory  (Ref.  2) 

(b)  Classical  Plate  Theory  (Ref.  4) 

(c)  Improved  Plate  Theory  (including  the  effects  of  shear  deformation  and 
rotatory  inertia).  (Ref.  4) 
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Tine  T  IN  SECONDS  »<10“ 

Comparison  of  I«P>7«  and  C«P*T«  with  Elasticity  Solution 


In  each  case,  an  exact  solution  (In  series  form)  was  found  for  the 


boundary  value  problem,  A  comparison  of  the  three  solutions  reveals 
essentially  equivalent  results  for  the  gross  motion  (deflection,  bending 
moment,  etc.).  In  this  connection,  see  Fig.  2.  This  implies  that  future 
calculations  for  practical  plate  (or  shell)  structures  subject  to  Laser 
irradiation  can  be  carried  out  within  the  framework  of  relatively  simple 
mathematical  models,  resulting  in  a  considerable  reduction  of  the  computa¬ 
tional  effort.  Fig.  2  shows  the  time  history  of  the  center  deflection  of 
an  irradiated  plate.  The  three  curves  correspond  to  the  three  different 
theories  employed  for  the  computation. 

In  references  1  and  3  we  consider  the  boundary  value  problem  of  the 
cJamped  circular  plate  subject  to  normal  Laser  irradiation  at  the  center  of 
the  plate.  An  exact  solution  for  the  dynamic  response  of  the  plate  was 
found  in  series  form.  This  model  includes  the  effects  of  flexure,  shear, 
transverse,  rotatory  and  radial  inertia  forces,  etc.  The  results  of  this 
study  served  as  the  theoretical  base  for  experiments  described  below. 

2.2  Experimental 

Thin  aluminum  and  steel  plates  in  the  shape  of  a  circle  (radius 
%  11.5  cm)  were  subjected  to  laser  irradiation  in  the  laboratory.  The 
plate  was  clamped  at  the  boundary  with  the  aid  of  a  heavy,  machined  ring 
fixture.  The  laser  beam  was  directed  to  intercept  the  plate  at  the  center, 
normal  to  the  plate  surface.  The  Laser  used  was  a  Holobeam  model  630-QNd 
glass  system.  This  Laser  produces  an  output  power  (in  the  Q  switched  mode) 
of  approximately  3  joules,  with  a  pulse  width  of  approximately  40x10"^  sec. 
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The  output  beam  is  about  1  cm  in  diameter  with  a  beam  divergence  of  2  m 
radians,  at  an  output  wave  length  of  1.06  p.  The  beam  deposits  radiant 
energy  onto  the  metal  plate,  resulting  in  heating.  Time  dependent  thermo- 
elastic  stresses  are  relieved  by  the  motion  of  the  plate. 

The  motion  of  the  plate  is  sensed  by  a  Bruel  and  Kjaer  model  MM004 
capacitive  transducer  and  associated  circuitry.  The  amplified  signal 
is  displayed  on  a  Tektronix  545B  oscilloscope  and  photographed.  The 
detector  was  mounted  at  a  separation  of  1  mm  from  the  rear  surface  (away 
from  the  laser) ,  at  the  center  of  the  plate. 

The  results  of  experiments  (plate  center  displacement  vs.  time)  are 
shown  in  Fig.  3.  The  theoretical  curve  is  also  shown  in  Fig.  3,  and  it 
can  be  concluded  that  the  mathematical  model  predicts  the  motion  with 
I  reasonable  accuracy. 

In  addition  to  the  above  experimental  work,  exploratory  work  was 
performed  to  move  the  plate  by  impulsive  means.  Lasers  were  used  to 
vaporize  a  very  thin  layer  of  the  plate  metal  surface  material.  This 
produces  a  plasma  which  is  heated  (by  the  Laser) ,  and  subsequently  pro¬ 
duces  a  shock  wave.  The  shock  wave  impinges  upon  the  plate  surface  and 
causes  the  plate  to  move  impulsively.  A  mathematical  model  for  this 
type  of  pressure  loading  has  been  considered  (see  reference  7) . 
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1 .  IVIKOI'UOTION 

It  is  well  known  that  lul cr - Betnou 1 1 i  beam  theory  and  its  two- 
dimensional  counterpart,  classical  j'late  theory,  neylect  the 
effects  of  rotatory  inertia  and  transverse  shear  deformations.  As 
a  consequence  of  this  omission,  results  obtained  hy  these  classi¬ 
cal  1  approx imat e '  theories  are  valid  only  for  the  case  of  waves 
whicli  are  lony  compared  to  the  radius  of  yyratii>n  of  the  beam 
cross- sect  ion  or  the  plate  thickne-is.  On  the  otlier  hand,  classi¬ 
cal  olastieity  theory,  while  st i Jl  limited  to  sufficiently  small 
deformat i ons ,  imposes  no  restriction  upon  the  wave  lenyth  ratio, 
init  solutions  for  I'cam  and  plate  dvnaiuics  prohlcnis  within  the  fra:iu  • 
work  tif  t  lu-  i.‘,\a,.  t  theorv  require  a  monui:K*ntaJ  ci'mpuf  at  i  ona  1  effort, 
if  such  solutions  can  he  ol>talned  at  all. 

In  the  case  of  hcams,  correct  ii'ii  terms  have  been  supplied 
hy  Ravleiyli  |1|  and  lim.'sluiko  (’.M-  Rayleiyh  introduced  the 
effect  of  rt>tatory  iiurtta  and  limoshenko  included  also  the  effect 
of  transverse  shear  deformation.  Ihus.  this  I'ne-d  i  mens  iona  I  thec'ry 
well  know/)  as  iiiimslienko  beam  ihei'r>'  occupies  a  |>osition  intermediate 
between  [  iiler-Bernoul I i  beam  theory  and  elasticity  theory  and 

"  J  Jl'  Si  I'Jt'J  1  i'35dO.  Jlh' 0  SM  Archives,  J  [lUSn  J13-J“" 
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is  I'Apoi'ti'J  to  (live  sat  i  sl'oel  ory  results  for  short  waves. 

Ill  tlie  ease  of  plates,  eipiations  of  motion  analogous  to 
Timoshenko's  heam  eipiations  have  been  yiveii  by  Miiuilin  |1|  aiul 
llflyaiiil  lri|,and  a  eorrespond i iiy  theory  of  plate  eipi i  1  i lir i urn  has 
been  yiven  by  lleneky  [bl  and  Keissner 

Sinee  these  improved  theories  and  eb  aentarv  tlieories 
are  eliaraeteri  eed  by  differiny.  dispersion  relations,  ue  can  expert 
that  eertain  wave  reflection  and  transmission  problems  associated 
with  piecewise  non- liomoyeneoiis  plates  and  beams  will  result  in  major 
differences  in  the  predictions  of  the  t  >  theories. 

The  present  investigation  is  divided  into  two  parts: 
jl)  Transmission  and  reflection  of  waves  in  piecewise 

non-liomoyeneous  beams.  In  this  case  we  consider  wave 
motion  in  two  bonded,  semi - i nf in i t e  beams  composed  of 
different  materials,  and  also  the  case  of  a  beam  of 
finite  lenyth  bonded  at  each  end  to  two  semi - i nf i ni t e 
lieains  composed  of  a  different  material. 

(d)  Transmission  and  reflection  of  waves  in  piecewise  non- 
homoyenoous  plates  the  construction  of  which  arc 
similar  to  part  (ll,  except  that  we  consider  wave  motion 
in  plates  instead  of  beams. 

1.  liQUATION  01-  BT.AMS 

2.1  Di'iUTtption  of  tho  Notii'ii 

We  shall  assume  the  following  displacement  components  for  a  beam: 

=  :  4Hx.t) 

ll  =  0  (  1  .  U 

y 

11^  -  w(x,t1 

Tlie  st  rai  n-d  i  splacement  relationship  can  I'C  written 
as  follows; 


X.N  “  3x 


1  1. 

‘•'.x:  =  3  V  *  J^j 


J.J  r>:t ionts  •Wii  iiarii  lion’s  Fninoiflo 

In  the  theory  of  small  motions  of  an  elastic  isotropic  contimium, 
the  kinetic  and  potential  energies  of  deformation  are 


T  =  \  PllV  ♦  U^,  *  iVjdx  dA 
=  /J'/^  P(z^  *  iv^)dx  dA 


=  3  4'4  ^^x'-'xx  *  VyV  * 


210  +210  +  2t  0  )dx  dA  (14) 

X3  :;3  XV  XV  vi  vz 


'  I  [\k'-  If  ^ 


where  a  dot  indicates  partial  differentiation  with  respect  to 
time.  Define  the  moment  of  inertia,  the  area  of  the  section,  the 
bending  moment  and  shearing  force  as 


1  =  dA 


A  =  dA 
A 


(1.6) 


J16 


H. 


M  =  /.  T  z  dA 
A  XX 


Q  =  ^A  ^xz 


Then  equations  (1.3)  and  (1.4)  reduce  to 


T  =  4  /*■  P(1  A  K^)dx 


(1.7) 

(1.8) 


The  variation  of  T  and  v  are 


5T  =  p(  I  (*i  6  ()'  +  A  w  6  k)  dx 


6V  =  /J  M  —  (6^1)  +  Q  |  .54,  +  i-  (6w)|J  dx 


(1.9) 

(l.iO) 


t'enote  the  work  done  by  surface  tractions  when  the  di.splaccments  arc 
varied  by  6W.  Then 


x=ll 


5W  =  (m6<^  ♦  q6w1dx  +  +  Q*5w] 


(1.11) 


where 


[t  zjr  dv 
V  ‘  xz  'h  ' 


V  t^zz'h 


M*  =  /,  T  z  dA  at  x=0  or  ? 

A  XX 

Q*  =  /,  T  dA  at  x=0  or  f 

A  X  z 

Hamilton's  principle  in  an  interval  of  t  ira<'  ti  to  t?  is 


/  (6T  -  5V  ♦  6W)dt  =  0 

t2 


I 


jr 


ii.i:) 


After  sulist  i  tilt  ion  of  ci|uations  (1.9),  (1.10)  and  (1.11)  into  (I.IJ) 
and  upon  application  of  inte!;ration  by  parts,  we  arc  led  to  the 
fol  lowing  ci|uation 


+  (-pT^  t  -  Qni)5i)i  I  d.x  dt 


/Vf.-PAW  . 

tl  o  >- 

t2r  -1X=II 

+  /  |(M*-M)iSi{)  * 


dt  =  0 


( 1  .  1  .A  I 


2.0  Equations  of  the  Motion 


From  equation  (l.l.i),  the  equations  of  motion  can  be  written  as 
follows : 


-PAW  *  g  *  q  =  0 

TV  II 

-p:<^  + 


(1.141 


The  approjiriate  boundary  conditions  are: 

11  At  each  end  of  the  beam,  x=0  and  x=t,  one  member 
of  each  of  the  pairs,  (M,(t')  and  (Q,wj  must  be 
specified. 

2)  At  the  surface  of  the  beam,  z  =  *  v*  one  member 
of  each  of  the  pairs  (ro,if)  and  ((|,wl  must  he 
speci f ied . 

To  insure  a  unique  solution  the  following  initial 
conditions  arc  required: 


w(x,il),  (f(x,0),  w(x,0)  and  if(x,(M  are  specified. 


case  to 


The  general  Hook's  law  reduces  for  the  present  isotropic 


i; 


V  l  I  T  . ) 

VV 


Ji: 


11.151 


Upon  sulist  i  tut  ion  of  equation.'^  (l.-l  into  (1.15),  we  obtain 


further  suiist  itut  ion  of  equations  (1.1(0  into  (1.(0,  and  usini; 
f  iiiioshenko '  s  eorreetion  factor  e,  tlie  beam  stress  displacement 
relations  become 


M 


li  I 


34. 

3x 


g  =  A  tl 


(i.ri 


where 

V  :  (T^,^,  t  t.JdA 


is  i>;nored. 

I'he  displacement  equations  of  motion  are  now  obt, lined  by 
sulist  i  tut  ion  of  (1.171  into  (l.M),  with  the  result 


■'  ‘  *  ‘I  =  PAW 

i;i  1^4  -  A  r.  +  4.^  +  m  =  pl0  .  (1.18) 


r(lijations  (1.18)  arc  the  well  known  Timoshenko  Beam  ei^iiations. 


Jl'.l 


I  \  >:on-hon.^:t:neo:i^^ 

J.4  F’.cv.ii-  F!:<x 

With  reference  to  equations  (1.17),  the  sum  of  kinetic  and  poten¬ 
tial  enerjjics  for  a  heara  of  lensth  xi  -  X2  is 

2  2  . 

T  *  V  =  /  ,.\W^)  .  i:i  leV  i  Jdx. 

(1.19) 


2  2 

/  ,.\W^)  i:i  Idx 

XI  L  ^ 


.•\fter  differentiation  of  (1.19)  with  respect  to  time  and  upon  appli¬ 
cation  of  integration  by  parts,  we  obtain 


S 


xz 


(TA-)  =  / 

.  ><> 

+  [Ml))  ♦  Qw] 

^2 


«(p.U.  -  g)  *  -  g  .  Q)]dx  . 


(i.:o) 


When  there  are  no  applied  loads,  q  =  m  =  n,  and  the  integrand  of  the 
right  hand  side  is  identically  tero.  Thus  the  equation  reduces  to 


-g  (T*V)  =  [M4>  Qw]^^ 


(i.;i) 


We  identify  the  energy  flux  as 


,J(x,tl  =  -  [Mij  Ow]  . 


(1.::) 


Upon  substitution  of  equations  (1.17)  into  (l.dd)  we  obtain 


.I(x 


.t)  =  -  El  If  i  *  ’"'''‘■(lie'  *  <>)'*■] 


(1  .-’31 


2.5  FeJut’tion  from  Timoshenko  Beiim  Hicoru  to  Fulcr'-Bemoui !  i 
Bern  Theory 

w 

If  we  let  0  =  -  Tp-  in  the  equation  (1.13),  we  are  led  to  the 

oX 

following  equation 


Tlio  dlsphu'emoiit  ei|uations  of  motion  and  moniont  and  shear  foree  are 
otiiained  as  follows: 


( 1 . do  I 


5.  KAVf  PROPAGATION  IN  PlfCfWISf  NON-IIOMOGPNPOUS  BPAMS 


3.1  ik'yicval  Sohtt'ion 


rquations  (1.18)  with  cj  =  ra  =  0  reduce  to 
*  *)  = 

[;r  -  c'AG  ♦  <{>^=  pi$ 

Assume 

K(x,t)  =  A 
4i(.x,t)  =  0 


(1.:-) 


(l.dS) 


where  1)  is  wave  frequency  and  C  is  phase  velocity.  Upon  suhstitutioi 
of  equations  (1.28)  into  (1.27)  we  obtain  a  pair  of  homoRcnoous, 
linear  algebraic  equations  in  A  and  U  whose  determinant,  set  equal 
to  :ero,  yields  the  velocity  equation 


(1-a)  0“  -  (Ce^  Cs^)  G^  +  =0  (1.21)) 


where 

(■?  =  li  (•’  =  1  =  ^  (-2 

e  r  ’  s  p  ’  '  isT^  ^  1  's 


Ironi  equation  (1.  'qi.  we  obtain  tlie  mode  veli>eities  as  fellows: 


f  ‘  *  C'  -  \  ((■■  +  l"  )  - 

e  s  e  s 

"  '  Jll-al 

C'  +  (“  +  I  (('*  ^  C*  I  ' 

e  s  e  s 


1 1  1  -a  )(“  ] 

e  s 

-  in 

-I  (  1  -a  i("  j 
e  s 


1  .  .'0  I 


TIu'  most  j;oiK'r:il  solution  of  is  tl-oroforo 


k=l 


k=l 


i-.(t 

ii^lt 


o  1 


1.  (t 


i-Ut 


I  1  .  k  11 


where  it  is  uiulerstood  tluit  the  frequeneies  must  I'e  the  s:.'me  for 
nil  waves  and  that  i:.-  will  hoeome  ima.cinary  as  a  ''  1.  A  plot  of 
phase  veU'vity  vs.  frequeney  i  .s  shown  in  1  iqure  1.1.  N'oti*  tli.it  of 
the  S  coefficients  appearinj;  in  et|uation  ll..'l1,  only  1  .ire  in¬ 
dependent.  After  suhstitution  of  equations  (l..il)  into  (l-d")  we 
find  that 


where 


A  =  P  P  i 
k  k  k 


«k  =  -  ''k'\‘  ’ 


C 

''k  =  "•—- 


k  =  1,: 

k  =  1,: 

k  =  1,: 


( 1  ..i.'i 


In  liuler-Hernoul  I  i  Beam  Theory,  eipiation  (1--('1  with 
ci=('  reduces  to 


i:i 


0..\K 


I  1  .  ■ 


SirK’C  l!;  is  iilw:iys  i nary ,  A:  must  1h‘  cu.n.il  to  zero  tv'  ij)Sure 
a  iuninUeJ  solution  for  \  *  -•  .  A  plot  of  plniso  volov.'ity  vs. 
frequonvy  is  shown  in  [•ij:uro  1.1. 


3..;  ivV;;’,'  :n  Tur  /-.-.i'w  ■  •’ 

Two  sonii-iiiTiiiite  beams  of  different  materials  are  bonded  at  ,\  =  (i, 
as  shown  in  I'iHure  l.J.  for  a  di  stiirlianee  eomini;  from  the  negative 
X  direetion.  let  w,  and  be  the  iiicomini;  waves,  w[  and  li  be 

the  retleeted  waves,  and  let  and  4^2  be  the  t *'ansm i 1 1 ed  waves. 
We  have 


Mlx,t)  =  A  e 
k=l  ‘ 


isilt  -  } 


>fi(.x.t|  =  r  I)  o 
k=l  * 


ifat  -  ) 

ki 


<  X  <  0 


Wi(x,t)  =  ;;  B. ,  e 


ifiit  *  jA-) 


*''(x,t)  =  i;  R  C 

k=l 


in(t  w  A-  ) 

kl 


WjCx,t)  =  L  A,  ,  0 

k=l 


iSUt  -  jA  ) 


«2(x.t)  =  T.  I),,  e 

k=l 


iS^(t  -  A-  ) 
kJ 


,  0  <  X  <  00 


t  1  .  .V.) 


A,  =  p,  .  n,  .  i 
kj  kj  kj 


B,  .  =  -  p  .  R  i 
I'.l  '  k  I  k  1 


C.  .  V/  . 

n  .  - _ 

'ki  i:(c!7AIPCi 

k.l  SJ 


k  =  1,’ 
.i  =  1,2 


where  the  subscript  j  =  1,2  refers  to  the  respective  domain 

of  the  lieain  tsee  I'ijjure  1.2),  and  the  subscript  k  =  1,2  refers  to  the 

mode  of  motion. 


l.F  -  hondeF',  .a— ,'f-j e/Yeftc  .~t\v~.s 


At  the  junction  .x=(>,  the  following  four  boundar)'  con¬ 
ditions  must  be  satisfied 


W;  (0 , 1 1  =  Wj  (0, t ) 

4>i  I0,t)  =  <>.-  10, t) 
Qi(0,t1  =  Qj{0.t) 
Mi (0,tl  =  Mj (O.tl 


where 


J  .1  \  ex  .1/ 


3$ . 


M,  =  h.  I . 

.1  J  .1 


3x 


and  where 


Oj  =  Wi+Wl 

+• 

=  ‘i’l  +  4>i 


W2  =  W2 
+ 

.‘'2  =02 


.1  =  1.- 


i  =  1,2 


(1  .”1 


Upon  substitution  of  the  eipiations  (l.,>hl  into  (l-o'^l.  we  olitain 
a  set  of  simultaneous,  linear,  aleebraic  e>("ations  (in  matrix  form): 


J  J(' 


Diicc  Ihii  1  >  incoming  waves  and  the  iiuiteri.il  pro¬ 

perties  arc  known,  tlio  ipiantities  Ki  i  ,  K.- j  .  l>i^  and  P,  ,  .ire  ol>- 
tained  by  the  a(>{>l  icat  ion  of  i!ramer's  Rule.  K)u‘n  a  ).  P.  i  must 
be  equal  to  zero  to  insure  a  bounded  solutivui  for  \ 

The  energy  flux  of  the  incoming  wave,  reflected  wave, 
and  transmitted  wave  is  defined  l\e  equation  eNcejU  that  on  1 

real  parts  of  are  used: 

Ji  =  -  I  (RcmP  (UoP  I  ♦  (Mh  UKoiy  I  1 

ji  =  -  ( ( uomO  ( KfO i  I  -  (IvqO  M 

.>*,  =  -  1  (RoMp  (Kc'.y  )  ♦  (RoQ*  1  (RoW*  I  1 


Tlu’  t  rarisini  ss  i  oil  ;iiul  rcf'cction  coot't'ic  iont  <  ,iro  thrretoro  obtiiiiH’i.i 
in  tlio  following  manner; 


•  hlO.i) 

■  Jlll'.tl 


1 1 .101 


v\horc  th*'  :iar  detiotos  the  time  averai;e  over  a  complete  perioJ. 

)ri>m  t)u‘  point  of  view  <^r  conservation  of  energy,  the  sum 
of  the  transmission  aiul  reflection  coefficients  must  he  equal  {o 
one,  i.e.  T  +  R  I.  Calculated  values  of  transmission  coefficie/it 
vs.  d  i  mens  i  on  I  ess  fro<.(Menc\'  are  [OtUted  in  I  i>;ures  1.3  throuy.h  l.h. 
In  the  i  u I er- Beriu'u 1 1  i  Beam  Ihei'rv,  the  d i  sp l act-ment 


equations  lu'come 

■x 

4* 

"i  =  'lit' 

iiUt  -  ) 

'  1 1 

. 

w  j  =  H  [  1  0 

iilt  +  A  (  i,.lt  -  ■' 

.  It,,  o  '•> 

(  1  .  11 

+ 

W;.  +  ,\  2  •>  e 

ii;|t  -  )  is;(t  -  y 

-  A,,  0  ‘'‘'■ 

1  0-  X-  • 

where  .. 

1  a''"’  /  1  1 

■  '-o  ■■ 

'  ■  hi  • 

J  A  , : 

At  tlio  junction,  x- 

0,  tlio  I'Ouiulnry  coiulitions 

"i (0, tl 

=  "■?  1 0 ,  t ) 

A"i (U.t) 

3w,  ((i,ti 

■■'Tx‘ 

Tx  ' 

(  1  .  IJ 

MitO.t) 

=  M:(0.t) 

Qitl'.tl 

=  QziO.tl 

wlierc 

32  K. 

M.  =  -  K. 
.1  .) 

*  i  •  j  '  1 .  - 

Q.  =  -  i;, 
^  i  1 

1.  ,  1  =  1.:: 

32  X 

yield  the  following 

sot  of  s  i mu  1  t.inootis ,  liiio.ir. 

11  l{;ol>r.i  ic 

oiiiniti<ins  (in  nintrix  form): 

I 


where 


Solving  the  equations  (1.43), 


)-i 

All  ■  (l^Y<FP^a(I+d-’) 

A  (1.14) 

”  _  :(ltYd^)(l*d) 

All  '  a{(i*Ya")’*2Yd(iHFT} 


Upon  substitution  of  equations  (1.44)  into  (1.41)  and  application 
of  equation  (I. 25),  we  obtain 


J*(0,t)  = 

1^1 1 

1  2 

.^1.1  ^ 

1 

.1^(0, t)  = 

|B,, 

1^ 

hill  ^ 

dtCO.t)  = 

1  Al  2 

1:2  >2 

C?2 

The  transmission  and  reflection  coefficients  are  calculated  to  he 


^  4(l^Yd=^)^(l^d)^Yd 
■I*((l,t)  {(l  +  YdM2  +  2Yd(l+d^) 

^  4Y^d^(l-d^)^^(l-Yd^)'‘ 
■ll((l,t)  {(l+Yd^)^<.2Yd(I+d^);^ 


I  1  .!(,) 


The  conservation  of  energy  check  is  satisfied: 


J 


('f  Fhtiti'  U'yi/.th  Bonded  at  Kaol-  Ftid 
dte-  Fi'^^irotu'd  ef  a  Pi  F^'eront  y*at,i‘:i'. 


A  hoam  of  infiiiito  loniith  is  borKlei!  at  each  end  to  two  semi  -  i  nf  i  Jii  t  c 
beams  eompo^ed  of  a  different  material  as  shown  in  1-i^iire  1.^. 

The  displacement  ecpiations  in  this  ease  are: 


^  '■<'  -  o) 

Wi  (>.  .t )  =  A  0  ^  '  kl  ' 

k=l 

.  -  ‘--(t  ^  '  ) 

.fi  (x.t  I  -  'Vi'-'  ^  I'l  ^ 

k=l  ^ 

h’l (X, t)  =  H  f  ^  '■kl  '' 

k=i  ''' 

|x,t  I  =  V  R  0  X  '  kl  ^ 

k=l  ''' 

H^-d) 

w,  (x.t)  -  A.  ,  0  '  k:  / 

k-i 

fiix,!)  =  ):  Ik  ,0  ^  'kj^ 


“  ’  t  X ,  t  ) 

=  V 

k=I 

e 

'  b: 

/  V- 

(x.tl 

k=l 

ii;| 

e 

"  ‘'k: 

W  3  (  X  .  t  ) 

=  T. 

ii2  1 
e 

t'b: 

1  X  ,  t  ) 

k=l 

-) 

==  5: 

k=  1 

K  0 

■’kA 

i  4L  1 

e 

(•  -  j  ■ 

0  ■-  X  ‘  f 


mm/A 

(i)  \ 

r2[  f2;  V 

Fijur'c  1.7  -  he.uin  of  Finite  LerijUi  FenJed  ic  Ttlc  i 
Beans  Coriposed  of  a  Different  .'-la  tenia  I 


The  relations  among  the  20  coefficients  in  equations  (1.181  are 
defined  in  a  manner  similar  to  that  of  equations  tl..>()1.  Sote 
that  C]  3  =  Cl  1  ,  C2  3  =  C2  1  ,  pi  3  =  Pi  1  .  P2  a  =  pa  1  ,  and  hi  3  =  hi  1  , 
ha  3  =  ha  1  for  this  specific  case. 

At  the  junction  points,  x  =  0  and  x  =  *• ,  the  follohini; 
boundary  conditions  must  be  satisfied; 


wi (0,t) 
❖  1(0, t] 
Qi  (O.t) 
Ml  (O.t) 
"'2  (S-.t) 
4i2(K,t) 
Q2(2,t) 
M2(f,t) 


=  “2  (O.t) 
=  $2 (0,t) 
=  Q2(0,t) 
=  M2  (0,t) 
=  «3((,t) 

=  ■liaCK.t) 
=  Q3(«,t) 
=  MaiS-.t) 


q .  =  .c; .  4*^  *  ■ 

.1  .1  j  j 

M.  =  i:.i.  M  . 

J  J  J  3a 


i  =  l.-*.n 

j  =  l,-\,3 


and  where 


4>i  =  4>i  3-  <l>i 


^■3  =  '(■a 


Upon  sut's  t  i  t  ut  ion  of  emiations  tl.tS)  into  (1.491,  wo  ont.iin  .i 
of  simultanoous,  llnonr,  alpcbraio  oiiuation.';  as  follows: 


1 

"u  '  i 


1  1  I  1 


k.  k.'  k  '  '  k.' 


,  i-i.  I 


''k'.'ki'ki  ' 


whore,  for  oxamplo 


fj  hj  G:  G,  ^ 

>  =  1-,-^  113  =  Gr  =  G,  • 


I’,  •  =  ,k-  p,  ■  =  ,  k  =  i.J  ,i  -  i,:,.4 

k  I  G,  .  -G  . 

k  1  s  I 


The  kpi.'int  it  ios  Kji,  K.- 1  .  Pi.’>  P.’.’.  Pi.’.  "is-  ■'"•G  >  aro 

obtainokl  by  tho  application  of  (iramor's  I’ailo. 

Tho  transmission  and  rofloction  oisot  fi  o  i  out  s  in  this 
aro  klofinoil  as: 


I 


■I,  .t ) 


T  : 


■I]  H'.t  I 


■  Is  (('.  t  ) 


I  !  . 


•  ii  (I'.n 


Klu'i'i'  Jjlx.tl  is  Jofinod  ill  ;i  iminiior  similar  to  oquations 
Two  luimorical  results  are  shoun  in  liguros  l.S  and  l.il. 

In  lai  1 er- Bernou 1 1 i  Beam  Iheory,  the  d i sn laeoment  equa- 
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The  correspond i ng  boundary  eonditions  are 
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and  uo  ol't^iiii  tl.o  followini;  s  inui  1 1  ;iin.‘cus  ,  linear,  :i  Ijjcl'ra  i  e 
oquat  i  ons  [  i  n  mat  ri  .\  t'orral  : 


''1  Pa 


TIu'  transmission  and  rot'lcetion  ooort'ic  i  ont  s  in  tliis  ease 
arc  dofinod  in  a  manner  similar  to  oi|i)ation  (l.I.s).  Tho  numorieal 
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1m  till.'  CMSo  ot'  I  u  I  or- Ik'iMUHi  1 1  i  Ik'.ini  llioory  tho  noJiir.i  is  dispersive, 
tVM'  modi  s  ot'  pro|);ij;at  ion  are  possible,  the  first  node  response  is 
alK.iv  s  a  travel  linvt  sinusoid,  and  voiile  the  seeond  node  response 
is  alwavs  an  attenuated  standiii);  wave.  A  study  of  Timoshenko 
lie.im  Theory  also  indicates  the  existence  of  tv>o  modes,  the  charact¬ 
eristics  of  v%hich  differ  from  that  of  foreitoiiii;  modes.  In  res¬ 
ponse  to  harmonic  excitation,  the  first  mode  response  |  lov^er  phase 
velocity)  v\  i  1  1  alv»ays  he  a  travcllini;  sinusoid,  while  the  second 
mode  response  Ihii’her  phvise  velocityl  will  he  in  the  form  of  a 
travelling  sinusviid  if  a-'  1  (high  frei|uenc  i  es !  or  an  attenuated 
standing  wave  if  a  ■  1  (low  f  rev)Uenc  i  es  I  .  It  can  he  sliown  that  the 
mean  eneru.v  flux  vanishes  in  the  case  of  a  standing  wave.  The 
corresjuniUing  phase  velocities  are  shown  in  I  igure  1.1. 

The  incoming  wave  is  assumed  to  he  sinusoidal.  With  reference  to 
figures  I..'  through  l.u,  we  note  that  limoshenko  Beam  Theory  trans¬ 
mission  coefficients  T  near  :cro  frev|uency  approach  those  of  lulcr- 
Bernoulli  Beam  Theory.  In  the  low  f rei|ucncy  range  we  have  aj  I, 
a;  >1,  and  therefore  transmitted  and  reflected  waves  have  Itravellingl 
sinusoidal  and  (standingl  attenuated  components.  In  the  high  fre- 
tjucncy  range  we  have  aj  <  1,  a,-'  <  1,  and  in  this  case  all  w.ive  compon¬ 
ents  arc  sinusoidal.  Targe  discrepancies  between  the  tvvo  theories 
occur  in  two  frci(ucncy  hands,  one  of  which  separates  the  high  from 
the  low  frei)ucncy  region,  and  the  other  of  which  is  the  high  fre¬ 
quency  region.  In  the  case  depicted  in  Tigures  1..T  and  l.h,  that 
region  is  characteriied  hy  m  ■'  I  and  a:  ''  1 ,  i.e.,  reflected  waves  are 
sinusoidal,  hut  transmitted  waves  are  sinusoidal  and  attenuated 


(  St  atki  i  !!>; )  .  rcsultiiii;  in  a  {'ror.'Uinvc-J  Iinvcriii'.;  rfu*  t  ran^rri  i  "  s  i  i  ti 

cvk-ft*  i  c  i  Lilt  .  In  r  ho  caso<  oxonqil  i  f  iod  !\v  !  ic^ros  1.5  aikl  i  . .  that 
roi'.ion  is  ^.bai'act  i-ri  o«.i  by  ai  '  1  aiul  a.  I,  and  in  thi''  case  t  lu' 
transmitted  »avt.‘S  are  sinusoidal  v^hile  the  ikfU-cted  uavi-s  are 
sinusoidal  and  attenuated.  In  the  hii;h  I'retiueney  !k->:ion  tr.in"^- 
riission  ooef'f ieients  ajffirckich  a  constant  \alue  uhivh  is  ](.)Utr 
th.in  those  of  [  u  1  er  -  Kernou  1 1  i  Beam  Iheory  as  the  t'retiueno  i  i-s 
apju'oaeli  int‘init\'.  It  is  also  interestiiiv;  tii  note  ( u  i  t  li  ikl'ereikc- 
to  !  li^ures  l.o  thiauiyh  !.(>'  that  in  tiio  region  uliieh  is  i  nt  er.ed :  at  e 
botvseen  the  U>vv  and  hijih  frc‘iuency  rei^ion  ti’.e  specific  variations 
of  the  transmission  coefficients  uitii  frequencies  are  strony])' 
de]''endenT  on  the  direction  of  travel  of  the  incannini;  i\ave.  ihis  is 
not  the  c.ise  for  tlie  1  uler-Bernoul  1  i  Beam  Theory  where  the  trans¬ 
mission  coefficient  is  i luk-pendent  of  t'requency  as  v\ell  as  t'f  the 
direct  Uni  of  travel  o\'  the  incoming  uave. 

4 .  0 

In  this  case  we  have  both  sinsuoidal  and  attenuated  v%a\k-s  in  all 
regions  of  the  beam  for  the  lou  freqiienc)’  region  (ai  '*  1.  a;  i 
and  only  sinusoidal  waves  in  all  regions  of  the  beam  for  the  biii'h 
freqra.ncy  region  (ai  ^  1,  a;  '  i),  |  or  the  particular  e\ari]Ue-  de- 

[licfcd  in  l  igures  and  I.!U  the  transit  jo/i  <kcar<  s.h.  n 

there  arc  simisoidal  waves  in  dv>main  'ai<  1  and  •' i  nuso  i  da  1  aiivi 

attenuated  waves  in  domain  (a.  ■  li  in  ligure  l._.  In  t:,-  Uu 

frequency  region  there  is  good  to  fair  agi\k::ient  between  the  tuo 
theories,  while  in  the  liigh  frequency  rcgiv'ii  there  appear'^  a 
difference  between  the  two  tlk‘orii‘<.  In  file  case  I'f  I  i  moshenfo 
Beam  llu'or)',  the  numbers  i>f  maxima  land  minimal  of  T  per  unit 
frequeric)’  approach  a  const  mt  value  as  the  J’lkquencies  in¬ 
crease,  while  in  the  case  of  1  uKr-Bernoul  1  i  Ik-am  llieory,  t  he\' 
decrease  as  tlie  frequencies  increase.  llu'se  diffelkuices  are  due 
to  the  difference  in  mode  |M*t>[H*Mies  and  p’.asi'  velocities  between 
the  two  tlieories. 
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PART  ri 

5.  BASIC  l-.QUATIONS  0[-  P1.ATHS 
5.1  Desjvivtior.  of  the  Motion 

We  shall  assume  the  following  displacement  ermponents  for  a  plate: 


Uy  =  :  ilJy(x.y,t5 

U.  =  w  (x.y,t)  .  C.  1  ) 


rise  strsiin-uisplacement  rclationshi|)  can  be  uritton  as 

fol lows; 


5.2  Enevay  Eonsidepation  and  Ham'  I  tan’s  Pp:  p!  a 


In  the  theory  of  small  motions  of  an  olsistic  isotropic  continuum, 
the  kinetic  and  potential  energies  of  dcform.it  ion  arc 


T  I  ■'  !  I  ;  HI  »ir  *11'  lil\-  Jv  cl;: 
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'  ^  1)  v> 


Jl  e  +Ji  o  )dx  dv  d: 
>■;:  y::  z%  ;:x 


a  j|,  I  j  \  XX  -X  tt  ?y  \  'v 

-  >  ->0 

*  'A-y  *  -y)*  '.X  {Z  *  -x)}'’'  ^ 

*  '^y(-y  *  ay)*  '^x(:'x  *  ■ 

whcro  a  dot  ind  ic-.itos  di  fforiMit  iat  ion  with  rospoct  to  time 
the  stress  resultants  are  defined  as  follows: 


M  =  1  :  t  d; 

X  XX 


M  ■=  d  1  r.  d: 

y  J„  y> 


:i4  •• 

^  h  -  i,,.L  ^  _  h 

*  f  " 

M  =  2  T  z  il  z  ,  :i  t  X  =  (I  o  r  X  =  ■- 

X  Jh  XX 
■> 

r  h 

M*  =  2  1  z  dz  ,  at  V  =  0  or  X  =  a 

1 

*  r 

M*  =  j  r  z  dz  .  it  \  -  0  or  X.  =  i 
xy  XV 


at  y  =  0  or  X  =  a 


at  X  =  (»  or  X  =  IZ 


at  y  =  0  or  x  =  a 


Hamilton's  principle  in  an  interval  of  time  ti  to  t?  is 


(.ST  -  ,SV  -  «K)  dt  =  0 


After  substitution  of  e(|uatioiis  (2.<i|.  (2. 'I,  and  IJ 
into  (2.9)  and  upon  application  of  intepr.it  ion  by  p.irts,  we  .ir 
led  to  the  following;  equation 


f  t  ?  /•  a  r  2.  f  /  ,  3  oM  9M  V 

L..  [  I  K-  iV  •  »r  •  -#  - 

•  (-  Tr  »,  •  -if  •  T?  ■  %  • 


3Q,  3Q,, 


+  I  -  phw  +  *  P  I6K>  dx  dy  dt  (2.10) 


,  I  1(M  -M  )6il;  +(M  -M  +(Q  -Q  )6w|  dv  dt 

ti  •' O  .X  X  X  XV'  XV'  '^V'  ^x  \x  ' 

...  JJ.Q 


f  ♦  ,  *  >'=^ 

*  1  dx  dt  =  0 

J  t  1  Jq  •  '  '  ■  ■  ■'  >  >  y  =  o 


5..V  Equations  oj'  the  Motion 


From  equation  (2.10),  the  ev|uations  of  motion  can  be  written  as 
fol lows: 

,3  ..  3M  DM 

.  Lh  ^  .  X  .  XV  ^  ^  ^  „ 

12  X  3x  Dv  ^x  X 


,  3  DM  DM 

.  Ph-^  V  ,  .  XV  _  ^  ^ 

1 2  Dy  Dx  y  y 

3Q  DQ 

-  "  3x^  *  S-y-^  "  " 


The  approf)riatc  boundary  conditions  arc 


I)  At  >{  =  0,  one  member  of  each  of  the  pairs 

(M  ,ii^  )  ,  and  (Q  ,w)  must  be  specified, 
xy  y  'x  ' 

J)  At  y=(),  v-a,  one  member  of  each  of  the  ria  i  rs  (M  ,v  K 
M  ,1*^  ),  and  (Q  ,w)  most  be  specified. 

XV  X  V  ‘ 


To  iiisoro  :i  iinii|Ue  solution,  the  followitii;  initial! 
J i t i ons  aro  rcqu i rod : 

,  ij.'^(x,y,0)  ,  C'^,(x,y,i>l,  '^^  ( .\ ,y , H) ,  K(.\,y 
and  w(x,y,!))  are  specified. 

The  i;eneral  Hooke's  law  reduces  for  tlie  present 
three  dimensional  case  to 
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Upon  sulist  i  tut  ion  of  equations  (d.ld'l  into  we  obtain 
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Since  the  second  terms  of  equations  (.’.III  are  neyliyible  in 


The  Jisphiccmcnt  ciiu.itions  of  motion  arc  now  ol'tainoil 
hy  substitution  of  (2.  In)  into  (J.ll),  with  tlio  result: 

V  ■"  3x3y  y*  ’  3>^  *  3x3y  J 

-  -  I'  V  S  =  y 

\ 3x  'xy  X  i:  ’x 

3x^)0* "  'y  (:ix3y  *  yx^O  ' 

-  (s-  *  f>  )*  %  "  1:  y- 

)f)*  K-c.hQ^j  *  ^  I’  =  .hw  .  i.'.n 


[■quations  iJ.l'l  were  olitained  ray  Miiullin  in  roferonee  ]  l|. 

3.1  .■■.':u' 

The  sura  of  Kinetic  anJ  potential  enerjjies  for  a  pl:ite  of  length 
X2  -  Xi  anU  width  >’2  -  yi.  with  reference  to  e(|uati>'ns  (.’..'1  aiul 
(:.41,  is 


*  V  =  *  y.  *  hw4 

-  jy  1  iL  l>-  '  y  I 

*  "x  •  "v  V/  *  \e('>y  *  if) 

*  '^xt.v  ‘  3'  )*  ‘V(fi  *  ’x)] 


After  d  1  f feront  i .It  i on  of  (J.ISl  uith  respect  to  t  inie  .u\d  upvMi 
.ippl  ic.it  ion  of  integration  l>y  parts,  ue  olUain 
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IVhoM  thoro  ;iri.'  no  appliod  loatls,  =  S,  '  I'  -  0,  t  lu'  tii'st  i  nt  opi'ai’.. 

of  the  rij^ht  hand  sido  is  i  (.lont  i  oa  1 1)'  zero.  When  tlu'  oi.l,dos  >  =  >i 
and  y  =  v;  of  tlio  plato  are  simply  supportod,  tho  third  i  nt  oyraiul 
is  idi-nt  ical  1\-  :oro.  I'hus  the  ociuation  rcdiioos  to 
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in  (.MUKition  [J.li)),  wo  aro  lo^I  ti'  the  tollowiiu',  t'nu.itioii 
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‘  •  J  x=o  J  y=o 

riuis,  the  equation  of  motion  can  !h'  written  as  follows: 
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X  ,  XV  V 

3P  3x3y  3y^ 


3h^/3‘w  S^iv'N  ^■'’x  ,- 

iy(,3P-*  3^>  -37  = 


If  \%c  Jclctc 


ph ^  ^ 

1 V  yy  yy 


•’S 

X  ,  \' 

ilfUl  -  ^ 


)■  t  “ 


wo  obtain 


3^M  3^1 ,  3’M  , 


3x^  “  3x3>'  3y 


r-  +  c  -  y"  +  f  =  i’l'W 
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SiiKA-  tlu'  fourth  aiul  fifth  terns  vani'^h  Incai-.-xc  of  the  eouatii'n 
(J.ll),  the  .i|’;Topr  i  at  v'  b,njnJar>‘  eojulitions  are  as  fol)<>us; 

il'  At  A  -  o,  x  ••  •»  one  r.erfner  of  each  of  the  pairs 
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must  he  sjK’Cified. 

( J :  \t  V  -  f>,  V  =  a.  one  nionhe/  of  each  of  the  jiairs 
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The.'  U  i  sp  1  acorai-nt  equ.it  ions  of  motion  and  !  ho  stre 
resultants  are  obtained  as  follows: 


b.  KAVI-;  I'ROI’ACATUIN  IN  P I liCblVI Si:  NON-IIOMOf.bNP.OUS  I'l.ATIiS 


('.1  -.'utioi! 

I'.quations  IJ.I'’)  with  -  P  =  0  reduce  to 


J.ll 
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-  (t)^ 


k  =  i,:,3 


IJ.SSI 


The  most  ^eiienil  solution  of  (J..i(i)  is  therefore 
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whore  it  is  nnderstooj  that  the  frec|uoncics  must  lie  the  same  for 
all  waves  and  that  can  become  imaginary.  A  plot  of  phase  velo¬ 
city  vs.  fre<|uency  is  shown  in  l  ipiirc  d.l. 

Soto  that  of  the  18  coefficients  appearinj;  in  ciiuation 
(2. Tib),  only  0  are  independent.  After  substitution  of  equations 
(J..i(i)  into  lJ..iO),  we  find  that 
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where 
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^  3  nt 


•k  IC^ 


,  k  -  l,d. 


>|3  =  " 


6.2  ^ave  Motion  in  Tuo  Bondcii  Stni^Inr'in:  tt-  /'/.iNi?  o.^td 
of  Different  Materials 


Two  semi- infinite  plates  of  different  materials  are  htindeJ  at 

x=0,  as  shown  in  Figure  2.2.  For  waves  coming  from  the  negative 

X  direction,  let  and  w,  be  the  incoming  waves,  v  v  , 

^xl  ^yl  1  ^  ^  ^ 

and  w,  be  the  reflected  waves,  and  let  'i»  ,,  ,  arid  w,  be  the 

1  ^x2  ^>'2 

transmitted  waves,  we  have 


ii  A*  c  kl  sin  —  >• 


S  ’  n- 

tl  =  ?.  B,  ,  o  kl  cos  -  >■ 

^,1  I'l  •* 


5  irat  -  -— ) 

w*(x,c,t)  =  X  B*  o  'kl  sin  --  V 

1  Kl  a 


3  _  iiat  * 

i.^,(x.y.t)  =  A^,  o  kl  sin  -  y 


3  _  is;(t  ♦  ~  1 

C.,.i(x.y.t)  =  B^l  e  kl  cos  --  y 


3  ifilt  ♦ 

wi(x,v,t)  =  t  nr,  c  kl  sin  --  y  I 

1  K=i  >  y 


(:.3i),'.) 

\  -.f  ^  X  "•  (1 


i:.3>lh) 


3  i Silt  -  -  ) 

i^'*i{x,y,tl  =  Z  A*,  o  kS!  sin  —  y  ] 


k=l 


3  illlt  -  .A  ) 

Kiyix.v.t)  =  ;;  B*  c  ^V.2  cos  —  y 

y-  k=i 


3  i:S(t  -  - 

w*(x,v,t)  =  y.  R*-,  o  kJ  sin  v 

k=l 


a 


(  J.3llc) 


where 


B,  .  =  p.  .  A.  .  i 
kj  kj  kj 


I'j 


■  k  I  k,i 


R*  .  =  q.  .  X*  .  i 
kj  ^kj  kj 


kj 


M.  ''r  ' 

k  I  k  I 


’’kj 


=  JlI  nn 

f:  a 


k  =  1  ,2 


■'kj 


_ ^ 

kj  K  (i. 


'’.J  =  -  :;jr 


k  s  1,2  .  q,.  =  0 

'•>1 


and  where  the  suhscript  j  =  1,2  refers  to  the  respective  domain  of 
the  plate  (sec  Figure  2.2).  The  first  subscript  k  =  1,2,.^  refers 
to  tl'e  mode  of  motion  and  the  second  subscript  refers  to  the  domain 
of  the  plate. 
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1 

> 

_ _ _ _ _ 

f 

) 

(5  ( 

;  ^2  ' 

f  • 

X«0 

Fi^jurc  2.2  -  Tuo  Bonded,  Semi-Infinite  Plates 


At  the  junction,  x=0,  the  following  six  boundary  con¬ 
ditions  must  be  satisfied: 

=  w,lO,y,t) 

=  M  ,  ( (I ,  V ,  t ) 
xy2  ■ 


il'jj ,  ( .  .V .  t ) 


Wj  ((),y,t) 

M  ,tO,v,t) 
X  1 


M  ,  1 0 ,  V  .  t  I 
xyl 


(2. .1(1) 


when',  for  o.\;iiiiplo 


and 


I  I’kl  I 


'■’kr'ki' 


Api 

An 


k=l 


7. 

k=l 


P,  ,  A 


kf  kl 


I'-  r,-. 


'ki''ki 


>r  i  "  r  '-  "  A 

n" 

a  '’k.i 

.  k=l,J,.A  ,  .i  =  l.J 

r  =  V  ,A- 

k  i  a  1.  . 

'’Ki 

,  k  =  1 , : ,  a  ,  ,i  =  1  . : 

k.i 

«ki  =  (t;  ‘‘kj 

+  1 

,  k^i,:,.A  ,  ;i  =  i,: 

Once  A* 1 ,  A2 1 ,  . 

A  3 1  and  i 

1  of  the  incoming  waves 

material  properties  are  known,  the  t|uantities  A,i,  Aj  i  ,  A31,  Aj.-., 
A22  and  A's  are  ol)tainod  iiy  tlio  application  of  Cramer's  Rule.  When 
Cii,  C?  1  or  C31  are  imaginary.  An,  A2 1  or  A31  must  he  ei|ual  to 
;oro,  respectively,  to  insure  a  hounded  solution  for  x  -> 

rile  energy  flux  of  the  incoming  wave,  reflected  wave, 
and  transmitted  wave  is  defined  hy  eipiation  iJ.dli,  except  that 
only  real  parts  of  |d..ii))  arc  used: 

I" 


d,=  -  I  (lKeM^|HRei;,^,IMReM^^,,l(Rej.^,j)MKe(J^,l(ReW|)|^  dv 
■'>’1 

/■ft  .  _  . -  -  •  -  ,  , 

.1;=  -  HReM'jHReiir|)MReM'^.,)(K<'>f^,j)*(ReQ^|)(ReW|)l^  dv 
Jy,  X  x  x,  .  . 

.Ij=  -  f  '  '  I  (ReM*  ,1  (Re.j,*  ,|+(ReM*^,  ,1  (ReiV.,l*IKeQ*  ,1  (Rcw’*;i  1^  dy  . 


Id.  l.’l 


transmission  and  reflection  coefficients  are  tliereforc 


ohtained  in  the  following  manner: 


R  = 


■Ij  10, t) 


d^'lO.t) 


■J,  (0,t) 


dj  (0,t) 


Id.. 13) 


where  the  bar  denotes  the  tine  average  over  a  complete  period. 

From  the  point  of  view  of  conservation  of  energy,  the 
sum  of  the  transmission  and  reflection  coefficients  must  be  ciiual 
to  one,  i.e.,  T  +  R  =  1.  Calculated  values  of  transmi  s.s  i  on  co¬ 
efficients  vs.  dimensionless  frequency 

./77h  , 


arc  plotted  in  Figures  2., 3  through  2.6. 

In  elementary  theory,  the  displacement  equations  liecome 

2  ifilt  -  , 

wt(x,v.t)  =  S  R.  ,  c  Hi  sin  y 

1  k=l 


w,  (.x,y,t)  =  j;  e 


k=l 


.  nr  , 
kl  sin  ■ —  V 


<  X  <  0 


where 


w,(x,y,t)  =  I  R,  T  e 
k=l 


in(t  -  /-) 

H2 


0  <  X  <  ‘O-- 

(2..1 ; 


'di  ■ 


^^4-  -  (T) 


At  the  junction,  x  =  0,  the  boumhiry  conditions 


awj  ((),>■, t)  :iw,  (d.y.t) 

'  ;h  '  ‘ 

M^j  ^  U'Ky.t ) 


>’ic!ci  the  t'ollou'injj  s  i  iiui  1  t;i  neons ,  line.ir  n  lyehr.i  i  e  etiiintioiis 


K=l,:,  .i  =  l,: 

The  eiierj;)'  fhix  of  the  incoming  wave,  reflecteil  iNOve. 
and  t-ansmitted  wave  is  defined  by  eoiiation  (J.dT)  except  that  only 
real  parts  of  id.-l-ll  are  tised: 


RT* 

sin  —  V 
a  ' 


Vxjtx.y.t)  =  V  e 


x-2 


k=l 


i;/_.,(x,y,t)  =  I  e 


k=l 


X  -  j? 

-  4T  - 

<*ns  — 


,  9,<x<°^ 


X  —  9. 

5  ^  iSi(t  -  ~) 

wtCx.y.t)  =  I  R^,  c  k3  sin  —  y 

k=l  ^ 


(:.48) 


The  relations  among  the  45  coefficients  in  equations  (2.48)  are 
defined  in  a  manner  similar  to  those  of  equation  (2.. 39).  Note 
that  Cj.^  =  Cj,j,  p^,.^  =  P|.j,  qj.3  =  q,.,,  ^  ,  f^,,  -  fj^j  and 

®k'  '  *^kl  '  1.2,3,  for  this  special  case. 

At  the  junction  points  x  =  0,  and  x  =  ^,  the  following 
boundary  conditions  must  be  satisfied; 


i('^l(0,y,t) 

'^^2(0,y,t) 

I'yiCO.y.t) 

= 

ij-y2(0>>'.t) 

wi (o,y,t) 

= 

W2 (0,y,t) 

Mxi(0,y,t) 

= 

Mx2(o.y.t) 

Mxyi(O.y.t) 

= 

Q,,(0,y,t) 

= 

Qx2(0"y.t) 

';'^3(2,y,t) 

= 

wj  (il,y,t) 

= 

W,(^,v,t) 

= 

Mx3(^,y,t) 

M  ,(l’,y,t) 
xy2 

= 

Mxy3(^,y,t) 

= 

Qx3(^.y.t) 

where 


uhoi'o  for  i  -- 


C.  .  i:  . 

K  «  •  k  I 


*kj  f.' 


nr 

'kj  "  <Yi  * '’o 


k- 1 


kj  :i  l',  .  '  k  j 


k=l,.',5 


«k  i  =  ‘‘k  i  '  ' 

k  I 


k=l  ..'.3 


.^1  =  c  .■  ,  -  V  (  .  , 

k  t,  k  kJ 

x,; 


k=) , J,3 


I  lio  t|mii\titios  A  j  j  »  A,’ 1  ,  Aji,  A^^,  A/^,  A^^,  A,.  .  A^;’,  ,  A13, 

aJ  3  ,  and  A33  aro  ol'taiiu'd  by  tlu-  appl  i  cat  i  on  iit'  Cramor's  ruU-. 

Tlio  transmission  and  rofloctii'H  cool'ricionts  in  this  oasr 
aro  dot' i  nod  as 

•itik.t) 


.1,  (d.t) 
d'lO.t  ) 


■  '1  ((),t  I 


whoro  Jjix.tl  is  doi'inod  in  a  mannor  similar  to  oipiation  (J.IPI. 
Nnmorioa)  rosnlts  aro  shown  in  I'iuuros  d.S  aiul  d.!i. 


//,  h  e 

wj  (n  ,y 

.t) 

=  K3a,y.t) 

3  Wj  {«  , 

y .  t ) 

3w3  (H  ,y  ,t) 

3x 

y.t) 

=  M  -(ll,v,t) 

xa  ■ 

'^xv:l 

3y  Jx=t 

r  3M  - 1 

=  ^xa  ^  Ilf  J  x=t 

and  wc  obtain  the  following  simultaneous,  linear,  algebraic 
equations; 


1  -1  1 

[  1  1 

1  1  1 

0  1 

Rii 

Rm 

'  ^^-1  ' 

1  ,1-  1 
k’ 

1  -  f--  1 

^k: 

<1  1 

It?  1 

rIi 

(  1 1 

[  -‘’kl' 

["'Pks  1 

1  mp|.,  1 

"  1 

Ri  ? 

Pi iRi  1 

f  "kll 

(mrik^,  1 

1  -ro'lj, ,  i 

0  ) 

R?2 

Pi ir! 1 

[  0  1 

f  ^k  ( 

(  , 

-i  1 

lit  2 

0 

1  0  1 

l(^  '“^kl 

^'k2 

-  1 
^•k3 

tin 

0 

1  0  1 

I'nPl^aS'l 

in.Pk/;  1 

-  ‘Vsl 

Rl  3 

11 

1  (1  ! 

'  "'^k2^k 

ll-n/k' 

1 

* 

R/3 

11 

- 

-J 

—  ^ 

r.‘> 

Numerical  values  of  transmission  coefficients  resulting  from  the 
elementary  theory  arc  also  shown  in  l  igiircs  J.8  ami  J .  H  for 
comparison  luirposes. 


CONCLUSION’S  AND  DISCUSSION  01'  RHSULTS  ll’AII'l’  III 


7.1  P,\}tes  of  Unboimicd  Lenjth 

In  the  ca.'^e  of  elcracntary  plate  theory  there  will  he  a  tr.iteliny 
wave  as  well  as  an  attenuated  standinj;  wave.  In  response  to  li.i  11:1011  i  i 
e.xeitation,  there  will  he  a  single  travel  ini;  sinusoid  if 


where  is  the  cut-off  freiiuency.  I'or  the  case  Os  fl'ere 

arc  no  travel  inn  waves  and  both  "waves"  will  lie  standing  and 
attenuated. 

A  study  of  improved  plate  theory  indicates  the  ewistence 

of  three  different  cut-off  freipiencies  0  •  s  c  ,  where 

S2,,,,  SI-,.,  il„.  satisfv  the  functions 
U.  .il.  JC 


Ik  = 


which  wc  can  obtain  from  eipiation  (J..ill.  In  the  freiiuenci' 
range  tliere  are  no  traveling  waves,  and  we  have  three 

standing,  attenuated  waves.  When  there  e.vists  a  single 

traveling  waive  and  two  standing  attenuated  waves.  Lor  ...j.-  :• 
wc  ha'e  two  traveling  waves  and  a  single  stamling  attenuated  wave. 
Wlicn  S.',j.  <Si  the  three  resulting  waves  are  of  the  traveling  varieti'. 
I’lots  of  the  phase  velocity  for  traveling  waves  .are  shown  in 
figure  .1.1.  A  dimensionless  plot  of  the  three  cut-off  frci|uencies 
vs.  plate  thickiies.'.  i.s  shown  in  I'igure  2.10.  Ilu'  cut-off  fre- 
ipiencies  12, and  become  unbounded  as  h  -D.  However. 


as  h  Ml.  fhese  will  corri'Siionil  to  the  vasi.'  of  id-assii'al  [ilate 
theory. 


7.2  Tuo  l^ondcil,  .'h-ni- Inj'initi' 

Materials 

The  incomin}!  waves  are  assumed  to  l>e  sinusoidal.  With  reference 
to  l-igures  2 .  .S  throuiili  2.(i  and  to  luKure  2.10,  we  note  that  trans¬ 
mission  coefficients  T  of  improved  plate  theory  .ippro.ich  those  of 
elementary  plate  theory  as  the  thickness  of  the  plate  approaches 
tero.  We  also  note  that  the  two  theories  are  in  aKroeraent  when 


In  the  hi^;h  frei|uency  ranj;e  Si,^.  £i,  T  for  improved  theory  approache 

a  const.'int  value  which  is  lower  than  the  T  for  element.iry  theory. 


:74 


/!.  Reit-fiaKK 

In  the  case  of  improved  theory,  the  transmission  co¬ 
efficients  for  each  of  the  three  modes  approach  the  same  constant 
value  for  sufficiently  high  frequencies,  and  this  value  is  inde¬ 
pendent  of  the  plate  thickness. 

An  incoming  wave  is  not  possible  for  In  the 

frequency  range  both  theories  predict  a  single,  travel¬ 

ing  wave,  and  the  transmission  coefficients  T  as  calculated  for 
both  theories  arc  in  good  agreement.  Differences  for  T  as  calculated 
by  the  two  theories  becomes  appreciable  for  <  £j. 

7.5  Pla'c  of  Finite  Lenoth  Bonded  to  Tuo  deni- Infinite  Fiatee 
Co>’i;'Osed  of  a  Vifferent  Motei'ial 

Kith  reference  to  figures  2.8  and  -.9,  we  note  that  there  is  good 
to  fair  agreement  between  the  two  theories  in  tlie  low  frei|uency 
range  <Q<n.^..  In  the  high  freciuency  range  there  can  tie  con¬ 
siderable  differences  between  two  theories. 

The  numbers  of  ma.xima  (and  minimal  of  T  (ht  unit  fre¬ 
quency  region  in  improved  theory  approaches  a  constant  value  as 
the  frequencies  increase,  while  in  elementary  theory  t  he.v  decrease 
as  the  frc<|ucncies  increase.  I'his  phenomenon  is  a  conse((uence  of 
the  difference  in  dispersion  relations  of  the  two  theories. 
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NOMENCEATllRi;  -  PART  11 

P.'in-iHpt  ior 

width  of  a  plate 
phase  velocity 

flexural  rijjidity  of  plate  (l>=l,h’/ 1 2  ( 1 -v^  ) ) 
Young's  Modulus 
exponential  function 


^'k.i 

M 

li 


coraporK'iit  s  of  strain  trnsor 

iiiojnlus  of  elasticity  in  shear  1  J (  1  + )  I 


^"(1 
e 

xy 

(i 

li  tliickness  of  a  plate 

i  -1 

,1  eneryy  flux  of  a  plate 

i  domain  of  a  jilate 

k  mode  of  motion 

I  lenyth  of  a  finite  plate  strip 

M  ,M  ,M  bendini;  and  twist  iny  moments  of  a  plate 
X  y  x\' 

M*,M*,M*  intensities  of  I'endiny  atid  twistiny  moment  at  I'oundar 
X  y  xy 

ra  ratio  of  Yoiiny's  Modulus  (ra  =  l.-'l  ;' 

1’  intensity  of  applied  surface  load 

Q^,Q^  shearing  forces  of  a  plate 

q*,Q*  intensities  of  shearing  force  at  houndar>' 

R  reflect  ioi)  coefficient 

S  ,S  intensities  of  applied  shearing  force  on  sui’face 

X  y  ‘ 

T  kinetic  energy 

T  transmission  coefficient 

t  time 

U  ,U  ,11  displacements  in  x,v  and  :  direction 

X  y  ;  ' 

V  potential  energy 

6K  virtual  work 

w  warping  displacement  of  plate  median  surl ace 

x,y,z  cartesian  coordinates 

Mindlin's  correction  factor 


Poi sson  '  s  rat  io 


di'Ms  i  r  y 

i.<^mj>onont  s  ot'  stress  tensor 

rotatit)fi  arit^ies  of  line  ori^einally  normal 
surf  aee 

f  re^|uene>’ 

the  kt!\  iiioJe  cut-off  frequency 
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1.  IMROinit  I  ION 

(he  two  prifie/pa)  oh  i  eet  i  res  of  the  present  report  ;ire:  111  to 
lo.uiealls'  d^telop  a  linear  theor>'  for  the  dyntimics  of  initi;ill>’ 
stressed  elistie  solids,  and  |Jl  to  present  ;i  vttriational  prinei- 
ple  witieh  will  seree  as  the  framework  for  the  systematie  develop¬ 
ment  of  approximate  theories  for  the  i nerement .a  1  motion  of  pre¬ 
stressed  rods,  beams,  plates,  shells,  ete.  Additionally,  the 
derivations  will  elearly  show  how  the  apptirent  meehanieal  jiroper- 
ties  of  the  solid  are  altered  hy  the  ])re'Stress.  To  aeeomplish 
these  ohieetives,  the  folhiwiny  assumptions  will  be  m;ide: 

111  The  solid  is  hypt'rel  ast  i  e ,  i.e..  it  possesses 
a  natur.il  state  in  which  all  stresses  and  strains 
vanish  s  imu  1 1  ani'oiis  1  \  .  and  its  neehtinieal  eonstitu- 
tion  is  eonipletel>  char. ict er i ted  by  a  strtiin  enei'B.v 
density  function  which  vanishes  ip  the  ntitural  st;ite. 
(J)  An  arhitrar>  [ire-stress  causes  the  initial  de¬ 
formation  of  the  solid  from  its  natural  state  to  the 
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;;o\  1.1'iu  n;.'.  the  liiKvir  thonry  i>t'  pro-st  rosseil  elastic  solids  is 
pros  on  t  id  .  I  h  i  s  rosult  prov  i  ilos  the  loiiioal  t'ramovtork  for  t  lu' 
..ti::,!'.  10  d  1 o  1  'pi’u'nt  op  soooial  apiM'i'x  i"iato  tlioorios  poi'  lu'o- 
st  I'issod  rods,  hiaiiis ,  plates,  shells,  etc.  A  resolution  of'  t  lu- 
I'i'sultinp  houndaip'  value  problem  I’l'r  the  ineromental  motion  op 
sol  ills  o  t'  homideil  extent  is  obtaim-d  by  tlie  application  oP  class¬ 
ical,  mat  heiiiat  i  ca  1  t  eclin  i  cpies  .  In  order  to  illustrate  how  the 
results  oi'  this  i  nvest  i  yyit  i  on  can  In-  applied,  two  specit'ic  cxaiaples 
are  presentcil:  (a  I  motion  oP  a  solid  subjected  to  an  initial  Ir.ilro- 
static  pressure,  and  (bl  the  motion  oP  an  initially  st  resscil  beaiii. 

Problems  relatin,r  to  the  behavior  of  pire-st  ressed  clast  iv 
soliils  are  op  considerable  interest  to  both  i-nyineers  and  physic¬ 
ists.  Inyinei'r,;  are  routinely  Paced  with  the  problem  of  analysin'.; 
,ind  des  i  >.;n  i  tu;  structtiri's  cont  .a  i  n  i  nc  [U'e-st  ri-ssed  members.  Ihis 
iiuariably  K'ails  tit  tin'  ipiest  ion  oP  st.ability  lor  instability'  oP 
the  initial  con  P  i  .yurat  i  on  .  Ihe  proper  way  ti'  assess  stability  is 
to  sul'iect  the  pro- st  r<.ssed  siWid  to  i  nc  roiiwiit  a  1  sorPace  t  mictions 
and  or  initial  cofulittons  atwl  to  obsec.c  wlu't.her  or  in>f  tlie  re- 
sultiir.;  incremental  motion  remains  close  to  the  initi.il  conPiyur.i- 
tion.  It'  it  does  then  tlu  initial  coni' i  y.iirat  i  I'n  i  -  stable  in  the 
sense  oP  l.ayranye  and  I  iapounov.  Ihe  closeness  rePi'rred  to  abo'.o 
must  be  iloPinei.1  ana  1  vt  i  c.i  I  I  y  in  terms  op  ,i  speeiPic  rwasiire  oP  the 
incremental  motion,  e.'C.,  the  s,piarid  norm  oP  the  i  nc  r(.'ment,i  I  dis 
pl.icement.  Moreover,  ip  the  (.'Nternal  I'lirces  which  induce  tin. 
pre-stress  are  conservat  i  \  e ,  then  the  initial  state  is  al'o  ^tahlc 
in  the  classic.-il  I  u  I  ei'- 1  ap.ranv.e  Sense,  i.e.,  the  (lotenlial  ener;;} 
ol'  the  s'stera  assumes  at  the  (.'(pi  i  1  i  br  i  ui:i  position  a  wx'.ak  i\'l;iti\e 
minimum  in  the  class  oP  virtual  d  i  sp  1  aei'incnt  s  sat  i  s  Py  i  ir.',  the  Pin- 
ematical  cuns  t  r.i  i  nt  s  .  \n  excellent  treatise  on  the  theoiu  ot' 
elastic  stability  In-  Knops  |||  appi'.irs  In  's .  I  luy.p.e's  I  my  c  I  ■ 'ped  i  .i 
')P  Physics.  He  i  n\  est  I  v.at  er;  both  I  hi'  I.  iapounov  criteria  and  i  he 
enerj;y  criteria  I'or  stability  .and  sbiiw:  when  t  lu'  twi'  .are  xuu  i  \  a  1  eiit 
atid  when  they  are  not.  .  liolotin  |  _’ |  has  also  i  n\ es  t  i  ya  I  id  ,1.ihilil\ 


!.ij  «.  tKier-iMin  -^nJ  1. 

from  the  dynamical  viewpoint  with  emphasis  on  the  study  of  non- 
conservative  external  forces.  In  addition  to  questions  of  stabil¬ 
ity,  the  study  of  pre- stressed  solids  has  been  pursued  because  of 
its  applications  to  the  field  of  crystal  physics.  The  material 
characteristics  of  solids  subject  to  initial  hydrostatic  pressure 
liave  been  used  to  test  various  atomic  model  potentials  for  crystals 
to  study  anharmonic  crystals,  and  to  study  various  thermal  and 
electrical  properties  of  crystals.  The  field  of  acoustics  has 
also  come  to  pla\-  an  important  role  in  the  present  subject  because 
most  of  the  experimental  techniques  used  to  study  the  material 
behavior  of  pre-stressed  solids  are  based  on  ultrasonic  methods. 

In  view  of  the  importance  of  the  applications  mentioned 
above,  a  great  deal  of  literature  on  the  subject  has  develoiied. 

Here  wc  will  cite  only  those  references  most  directly  relevant  to 
the  present  treatment.  A  more  extensive  list  of  references  may 
be  found  in  [1]  and  also  in  the  treatise  on  non-linear  field 
theories  of  mechanics  by  Truesdell  and  Noll  (.1].  As  pointed  out 
in  [.1],  Cauch)'  successfully  derived  ;i  line.ir  theorv  of  incremental 
motion  superimposed  on  ;i  pre-stressed  el.astic  bod\- .  fauch\'s  eq¬ 
uations  arc  given  l)y  fodhunter  ;ind  Pearson  l-l).  Sa  int -tenant  |1| 
attempted  to  duplicate  Ctiuchy's  results  using  an  energy  principle. 
However,  his  results  were  not  valid  because  he  did  not  retain  the 
quadratic  terms  in  his  expansion  of  the  strain  energy.  More  re¬ 
cent  derivations  of  the  linear  theory  for  incremental  motion 
superimposed  on  finitely  deformed  elastic  bodies  were  given  b>- 
Toupin  and  Bernstein  ( .S  |  ,  Truesdell  and  Noll  |  .S  |  ,  Hreen,  Kivlin 
and  Shield  |J8|,  (Ireen  and  3ern.i  |2;)|,  Knops  |1|,  l.r ingen  and 
Suhubi  |()|,  Of  these,  the  one  which  is  probablx’  closest  in  spirit 
to  the  present  vari.it  ional  ajrproach  is  the  derivation  given  b\' 

Knops  biased  on  the  invariance  of  the  rate  of  work  equations  under 
arbitrary  rigid  body  motions.  The  results  of  1 28 1  and  |  2P  |  are 
identical  to  the  ones  in  Section  2  of  the  present  i  nvest  i  g.at  i  on , 
but  they  .are  obtained  in  .an  entirely  different  manner.  Applic.ation 


of  these  results  to  particular  cases  of  initial Iv  stressed  solids 
can  be  found  in  1 30 1  throuj;h  l-iS).  We  also  mention  the  contri¬ 
butions  of  liiot  which  are  suiiiniar i zed  in  13(>1.  In  this  latter 
treatment  of  the  initially  stressed  solid,  v.iriahles  are  tailored 
to  tile  sjiecific  asymmetry  of  the  physics  and  jjeoiiietry  of  the  pro¬ 
blems  at  hand.  In  many  cases,  such  variables  are  not  tensors, 
by  sacrificing  the  requirements  of  invariance,  this  work  does  not 
seem  to  fit  into  the  main  stream  of  modern  continuum  mechanics. 

2.  AN.M.V.SIS  ,\.\1)  V.ARl.VnONAl.  1-0RMIII.A|- 1  ON 

Throughout  the  forthcoming  analysis,  standard  fartesian  tensor 
notation  will  be  used,  with  all  particle  coordinates  and  tensor 
components  referred  to  the  same,  fi.xed,  rectangular  tairtcsian 
coordinate  system.  Three  different  configurations  of  the  hi’por- 
elastic,  solid  body  must  be  distinguished: 

(1]  In  its  natural  state,  the  body  occupies  the 
region  V,^  bounded  by  the  surface  .S  ^  and  the  co¬ 
ordinates  of  a  ty|iical  niateri.i!  particle  I'  are 

<.lcnotc<J  bv  a..  I'hc  stress,  strain,  anil  strain 
•  1 

energy  dcnsit\'  all  vanish  identically  tiiroughout 

(2)  In  its  initial  (static  equilibrium)  state,  tlie 

body  occiqiios  the  region  bounded  by  the  surface 

S  and  the  coordinates  of  I'  are  denoted  bv  1  he 

o  •  I 

fauchy  stress  tensor,  l.agrangian  strain  tensor,  and 

strain  energy  density  are  denoted,  respect i ve 1 v ,  bv 

1^.’.,  i:‘?.,  and  K*. 

I  j  ’  I  j  ’  o 

(.3)  In  the  current  (or  fin.il)  state,  the  body  occupies 
the  region  V|.-  bounded  liy  tlie  surface  .S|,  .ind  the  co¬ 
ord  in.itcs  of  P  are  denoted  by  x*| .  Ibe  tlauchy  stress 
tensor,  l.agrangian  str.,in  tensor,  .and  strain  energy 
density  are  denoted,  :  espect  i vc I v ,  bv  t!  .  ,  f*.  .  ,  .and  IV,*. 
Ihe  initial  deformation  carries  the  material  particle  P  from  its 


I 


Uuntiim  i'.j  ( ,1  .  ,  :i ,  ,  ;i  .  )  in  the  initiirnl  f  i  sMunt  i  on  ii‘  tin-  li 

O  ... 

t  i  on  I'  l\  ,\  ,\.l  in  t  Ik‘  initinl  coni  i  I'ur.it  iini .  lln-  i  nc  I'li;:!. 

"  I  I 

lie  I  011,1, it  1  on  t  iu-n  enri'ies  it  to  its  linni  loention  ij  i\  ,i 

I'Ik'  eoiiiponciit  s  oI'  the  i  ne  reinent  n  I  J  i  sp  1  .lei-iiient  veelor  i’ ^l'|  si 

he  lieliot  I'll  I"  11 .  ,  i  .  e  .  , 

I 

I  o 

1 1  ,  =  , 

I  I  I 

I'.iti,  ii  pent  to  the  natur.il  eon  t' i  s’.n  rn  t  i  on  ,  t  Ik-  initial  aiul  ei 
I  on!  1:111. it  ions  are  ileseriheil  h\  the  oni--to-oiK-  iiiapp  i  nr.-' 


o  o 

\  .  .\  .  1  a  .  ,  a  ,  a  i 

1  I  ■  • 

I  1 

.  V ,  I  a  ■  ,  a  ,a  ,  t  ) 

1  1  ■  ' 


lon\ersei\,  ulth  respeet  to  the  initial  eon  f  i  ;;nra  t  i  on  . 

,  "  "  i-' , 

,1  ,  :  a  ,  X  1  ,  .\  .  .  I 

I  1 

aiul  t  hei'i.  tore  iJ.Jhi  ean  also  he  uritti-n  a-. 


I  I  o  o  o 

.X  .  -s  \  .  I  X  1  ,  X  ,  X  .  ,  t  I 
I  I 


Sinee  the  results  ot'  any  men  sn  reiiieii  t  s  perl'oriiieil  on  a  pri  -  t  r- 
hoi.l\-  are  ohtaineil  in  terms  ot'  the  initial  eoori.1  i  nat  es ,  Ue  -  h, 
11 1 1  i  mat  e  I  express  all  results  in  tei-ins  of  tlu-si.-  eooril  i  nat  i,  s 
us  i  It!.'.  I  J .  .-i )  aiul  I  J .  11. 

the  lay  rani;  1  an  strain  tensor,  ('.iiielix  st  re  ■  tin  or. 
eoiist  1 1  lit  i  \  e  relations  are  Jel'iiu-il  h>  the  t'a-iiiliar  ri.l.ition 
for  example  fhapter  ot'  1  ^  |  : 


le.i  - 

■nt  .1  I 


ii  ri  nt 


mil 

I  t. 


In  1  .  h  I  ,  I 


1 the  -.yminetrie,  I' i  o  1 ^  |  j  fehho  t't'  lie 


n  o 


i'  till,'  '.tr.iin  oiK'rv'.x  dt-niilv'  l'uiu'li"n.  i..-.,  I  !■■ 

[HT  unit  voluiik  ot'  t  iiu  n.ituiMl  In  .1  ' 

I ' 

.1,  -  Ji.'(  I  '.'i  I  i  ■-  tiu'  i  in  Ji-rir-', ;  nui:  •')  tin  iri 

‘  '  I 

lljH'n  sul'-t  i  t  ii!  ii>ii  nC  :  .  I  '  ililii  i-.r  !  i  n.i  I 

tutul  struin  uan  hr  u  r  i  1 1 1'li  .r.  t  hr  sui:i  ,it'  t  hr  iiMti.il  t 
an  i  nr  rrtirnt  a  1  st  ra  i  n  , 


i  -I  t  hr  initial  tt  na  i  n  .itui , 


s  t  hr 

i  nr  rriiirnf  a  1 

1  ''t  »M  i  f)  rt  t'rl  Trd  t  i>  t 

!u-  i  n  i  t  i  .1  1  t  .it 

s  iiSrJ 

t  r  lUauit  I,- 

p:i  ft  i  ;i  I  d  i  r!\  i  .it  i 

1  oil  u  i  f  ii  v'l 

nit  ial 

ruofii  i  nat  i 

.  i  .  f  .  .  -n,  '  •  \ 

)  K 

IVr  now  as-.iiiiir  that  t  lu  <tIMIii  rnrir.'.  ilrir  i  I  \ 


turrtiu'i'  kith  i  t -t  pai'fiaf  Jri'i  \  .it  i r  •  uj'  Ir  .ukI  inrIuJin 

thil'd  rr'ilrr  in  strain.  is\i'-t  aiul  arr  r.inl  inurii  ,  ,it  I  ' 

I  I 

statrl.  In  this  rase  uc-  ran  usr  larlrr's  i  arrrrri  tr  :i|'| 
in  the  nr  i  r.hhdiin "  si  nf  t  hr  initial  -t.iti-,  |>y  a  I'nnrtirn 
ilinulratir  in  t  lu  i  nr  rriiirn  I  a  I  -train,  i.-,, 

l\  ‘  I  I  '  .  I 
I  j 
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or  ill  \  \  t'u  of  ( .! .  ‘J.i  I  , 


In  view  of  this  result,  the  total  strain  energy  of  the  (hh(>’  is 
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is  the  strain  energy  associated  with  the  initial  deformation,  and 
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is  the  incremental  strain  energy.  To  obtain  these  results  we  made 
use  of  the  relation  d\'  =.l  d\  to  transform  the  volume  integrals  from 
the  natural  to  the  initial  configuration. 
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Ke  shall  now  proceed  to  derive  tlie  linear  equations  of 
incremental  motion  for  the  initially  stressed  solid  using 
Hamilton's  i'rinciple.  As  sliown  in  I”!,  Hamilton's  I'rincijile  for  a 
liypore  last  i  c  solid  takes  tlie  form 
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the  current  state  of  the  hod>'.  Similarly,  I'p  and  H|;  are  the  total 
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In  the  derivation  of  the  second  eiiiiation  in  (d.ddal  and  (d.ddh)  we 
made  use  of  iiauss'  Tlteorem  to  transform  from  volume  to  surface  in¬ 
tegrals  and  we  also  ira]iosed  tlie  eiiu  i  1  i  hri  um  conditions  on  the  initial 
conf  igurat ion ,  i . e. , 
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Tho  use  oF  (2.J2b),  and  (d.^rj)  in  (’..>"1  >'iclds  the  dis- 

placomont  equations  of  motion  and  associated  houndar>'  conditions 
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where  we  have  util  iced  (2.-4). 

These  equations  have  the  same  form  as  (1.T4)  and  (1.A5), 
page  -l(i  of  |,1|.  In  the  case  of  a  homogeneous  body  subjected  to  a 
constant  pre- stress,  (d.40)  becomes 
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fhis  form  is  identical  to  equation  (iii),  page  84  of  (4J  and  is 
attributed  to  Oaucliy  (,1829).  The  system  of  equations  and  boundary 
conditions  (2..T8)  also  appear  in:  jl]  (]>.  180,  (28.17),  (2.8.18), 
(2.8.19));  [81  (p.  247,  (<i8.9));  and  1(>|  (p.  288,  (4.2.80),  also 
sec  p.  2.84,  (4.2.8.11  compared  to  our  (2.4011.  The  present  deri¬ 
vation  differs  from  these  otliers  in  that  it  is  based  on  a  varia¬ 
tional  principle  (2.87)  which,  with  suitable  restrictions  on  the 
class  of  allowable  displacements,  can  be  used  to  derive  special 
appro.ximate  theories  for  pro-stressed  rods,  beams,  plates,  and 
she  I  Is. 
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cremental  deformation  wliich  is  compatible  with  the  boundar>'  con¬ 
ditions.  Knops  |1|  also  obtains  (2.11)  by  re(|iiirinK  that  the 
potential  energy  be  a  weak  relative  niiniimim  at  the  initial  ecpiili- 
briiim  state,  lie  refers  to  this  as  lladamard  infinitesimal  stability 
since  (2.11)  necessari  !>■  implies  th.it  lladaniard's  inecjualit>' 
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furtlier  shows  that  for  conservative  loads,  (2.  Ill  is  necess.i|-\  Init 
not  sufficient  to  guarantee  fiapounov  stability.  for  examiile.  in 
tlie  case  where  the  initial  eipii  librium  state  satisfies  displacement 
or  mi.xed  bound.iry  conditions,  the  inequality 
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for  all  symmetric  tensors  .  .ind  some  positive  constant  d  is 
sufficient  to  guarantee  l.iapounov  stability.  If  the  initial  state 
is  maintained  s(.)lel>'  by  traction  boundary  conditions,  then  in 
addition  to  (2.1(1 1,  the  restriction 
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is  necessary  to  prove  l.iapounov  st,ibilit\ 
these  cases  is 
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If  tlie  criterion  for  infinitesimal  stabilitv  (2.  Ill  is 
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satisfied  then  one  can  easily  show  tli.it  (2..SH)  possesses  it  most 
one  solution  if  the  initial  displacement  and  velocitv  are  specit'ied 
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t  h  rouslu'iit  \' ^  1 .1 1  .  I'urtlu'rmore,  the  natural  frequencies  assuci- 

ateil  iiith  free  liarmonic  incremental  motions  (f.=(),  T.=H,  1  ,=0I  of 
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tlie  form  u.=ll,i\)cos  t  are  real.  I'o  prove  this,  multiply  (r.-^Sa) 

t'\  II  .  integrate  over  \  aiul  impose  the  homoyeneous  houiuiar\  data 
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I  luis ,  ii.  In  tlie  case  of  strict  inequal  it>.  tl-O,  (referred  to 
as  superst ah i  1  i t >■  in  w'  further  ciuicliule  that  no  natural 

f requeue)'  can  vanish. 

An  energy  continuitv  eipiation  can  he  derived  for  the 
incremental  motion  hv  imposini;  the  conservation  of  eneiyu)'  law  in 
the  following  form: 
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and  e*  is  the  specific  energy  of  the  hody.  Suhstitution  of  (d.dSl) 
(J..tJl,  and  lJ..'ill  into  (J..s(M  yields 
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The  equilihrium  and  hoiindar)  conditions  for  the  initial  state 
tojdther  with  |J..SScl  may  he  used  to  reduce  the  right  hand  side  of 
I -’..Adi  to 
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s.  is  oalloil  the  oneri'v  flux  vector  siiiee  it  is  seen 
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the  eiierj;\  I'lou  per  unit  area,  per  unit  t  iiiii-  out  of  t  lu' 
The  application  of  liauss'  I'heoreiii  to  t  !u'  surface 
then  yielils  the  enerr.y  ia)ntiiniity  ecpiation. 
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Srateil  in  woials,  tlie  local  rate  of  increase  ol'  enerr.y  is  equal  to 
the  rate  at  which  the  hoih'  forci'  and  initial  stress  do  work  minus 
the  local  rate  of  et'flux  of  i-nerr.x  . 
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One  of  the  central  prohleias  arisin.r  in  enr, i neer i its;  applications 
of  the  thec'rs  dexeloped  in  Secti.ni  d  is  the  t'orced,  i  nc  r<.-rient  a  1 
motion  of  a  hininded,  initially  stressed  si'lid  with  t  iru-  ds'peiulcnt  , 
incremental  forciiyu  funct  ii>ns  and  Innindary  conditions.  Uc  shall 
resolve  this  pri'hlem  in  the  manner  <'l'  reiAu'enci.'  |  |  usins; 

classical  mathematical  techniques. 

Kith  reference  to  the  de\  e  lopment  in  Section  .  ,i  well 
posed  fs'rced  motion  problem  can  be  stated  as  t'l'llows. 

Iquation  of  incremental  iiu'tion: 
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It  can  bo  shown  that  the  solution  ot'  tho  probloiii  as  olia  raot  i.!’ i  ■  od 
by  oi|uation.s  (.i.ll  is  unii|uo  if  Id.  Ill  is  satisfiod  (soo  |  .i  1  .  |i. 
diil'  )  . 

To  olifain  a  solution  of  t.i.li,  wo  jirooood  in  a  mannor 
similar  to  Id'],  i.o.,  wo  assuiiio  a  solution  of  l.i.ll  in  tho  t'orm 
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where  6  is  the  Kronecker  delta  svmbol.  liquation  (5.51  also 
mn  (in) 

implies  that  the  eii’enfunctions  U.  have  been  suitably  normalited. 

In  view  of  (5. Id),  (5.1e),  {5.5c),  (5..5d),  (5.)cl,  .and 
{5.4d)  the  assumed  solution  (S.d)  satisfies  the  non-homoqeneou.s 
boundary  conditions  (5. Id)  and  (5.1e).  To  determine  the  scalar 
functions  4jjj(.t)  we  now  substitute  (5.:)  into  (5.1a)  and  (5.1b|. 

If,  in  addition,  we  apply  (5.5a),  (5.5b),  (5.4al,  and  (5.4bi,  we 
obtain 
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m=l 
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Similarly,  if  we  set  t  =  0  in  (5.:)  and  substitute  tlie  resulting 
equation  into  (5. If)  we  readily  obtain 


I  u["’^  (.x)*q^(o)  =  u.*°'(xl  -  u  '^’(x,oi 
m=l 


(  5 . bb  I 


Z  u['"^(x)*q^(o)  =  u/°*(x)  -  u.'^’(.x,o) 
m=l 


( 5. be  ) 


Now  we  multiply  (5.6a),  (5.bb),  and  (5.6c)  by  u|”*,  and 

respectively,  and  integrate  the  resulting  equations  over 
the  volume  V^.  Upon  application  of  (5.5)  we  obtain 
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VnIuTl'  tlK'  supiM'So  r  i  I't  a  n-t'ors  to  llu-  natural  stato.  Ii|uatioii 
^  I .  t' I  vorit'ii'S  our  assertion  that  onl>'  terms  up  to  the  t  li  i  lal  order 
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eorreet  ion  term  for  A'.*.,  .  Subst  i  Tut  i  in;  i  l.ldl  into  i  1  .  I  and 
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Hero 

lie  are  adopt  ini;  the  notation  of  loupin  and  Hernstein  |  .s  |  .  Substitu¬ 
tion  of  (l.lal  and  (l.“al  into  (l.lf)  yields  the  relation 

p  c,  .  =  il-di  f''-  I-  A\‘.,,  -  .  '.V.'.,,  I 
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and  since  we  are 
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lhi<  -^aiiK.'  rosult  is  also  obtaiiioil  I'.o  Toui'in  anil  Itornst  o  i  n  |  a  |  .  a-i 
uoll  as  In  riuii'stoii  aiul  Hru;,’,;K'r  In  an  i-arlicr  papor.  Hi  rob 

|lol.  ilcr  i  toil  a  rostriotoil  t'onii  of  il.JOi  basod  on  t  lu-  assumption 
that  h'  is  a  ipiadrat  i  o  t'linot  ion  ok  the  lailorian  strains,  rather 
than  the  1  ,i,i;raipp  i  an  strains.  In  terms  of  strain  invariants. 

Hireh  assumed  that 
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Hh'  substitution  ot"  (1.2.!)  into  (4.^11  yields,  after  neg  I  ect  i  nj; 
terms  higher  than  the  third  order. 
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On  the  other  hand,  substitution  of  (l.S)  into  (4.4)  yields  the 
resu  1 1 
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Birch,  lias  therefore  assumed,  instead  of  three  independent  third 
order  constants,  the  values 
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•'ll 


(4.:4) 


ihese  values,  when  substituted  into  (l.dlM  yield  Hrich's  results. 
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Ihese  formulae  are  very  appealing  since  the>'  contain  only  the  first 
order  elasticities  4'*  and  |i''  and  the>'  also  predict  that  the  wave 
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spc'i'ils  i  IK' ri'iisc  Kith  lirt'ssuro  as  oiu-  iiiipht  expect.  lloKever, 
recent  1>  Soca  and  Andersiin  |  1  .A  |  ,  have  shouii  experiiiienl  a  1  1  that 
i'i>r  i  sot  ro|i  i  c  ,  po  1  >c  rvst  a  1  1  i  ne  dnO  speeiiiieiis,  the  sliear  wate 
sfieed  decreases  Kith  incrt-asini;  pressure.  this  result  can  not  be 
predicted  usiin;  Birch's  f'orriiulae  I  l.J.'i).  lloKever,  the  i;eneral 
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inan>  authors  still  use  tlie  earlier  results. 
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